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Large Nongaussianity in Axion Inflation 
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The inflationary paradigm has enjoyed phenomenological success, however, a compelling particle 
physics realization is still lacking. The key obstruction is that the requirement of a suitably flat 
scalar potential is sensitive to Ultra- Violet (UV) physics. Axions are among the best-motivated 
inflaton candidates, since the flatness of their potential is naturally protected by a shift symmetry. 
We re-consider the cosmological perturbations in axion inflation, consistently accounting for the 
coupling to gauge fields (f)FF, which is generically present in these models. This coupling leads to 
production of gauge quanta, which provide a new source of inflaton fluctuations, 8<j>. For an axion 
decay constant ^ 10 -2 M p , this effect typically dominates over the standard fluctuations from the 
vacuum, and saturates the current observational bounds on nongaussianity of the CMB anisotropics. 
Since sub-Planckian values of the decay constant are typical for concrete realizations that admit a 
UV completion (such as N-flation and axion monodromy), we conclude that large nongaussianity is 
easily obtained in very minimal and natural realizations of inflation. 

PACS numbers: 11.25. Wx, 98.80.Cq 



I. INTRODUCTION 

Primordial inflation is the dominant paradigm in cur- 
rent cosmology since (i) it resolves the conceptual diffi- 
culties of the standard big bang model, and (ii) it pre- 
dicts primordial perturbations with properties in excel- 
lent agreement with those that characterize the Cos- 
mic Microwave Background (CMB) anisotropics. De- 
spite these successes, there is still no compelling par- 
ticle physics model of inflation; the key obstacle being 
the requirement of a sufficiently flat scalar potential, 
V (</)). Even generic Planck-suppressed corrections may 
yield unacceptably large contributions to the slow roll 

parameters e = -f- (jr) , f] = M^X^, thus spoiling in- 
flation (prime denotes derivative with respect to <fi, while 
M p = 2.4 • 10 18 GeV is the reduced Planck mass). One 
of the simplest solutions to this problem is to assume 
that the inflaton <i is a Pseudo Nambu Goldstone Bo- 
son (PNGB) In this case the inflaton enjoys a 
shift symmetry <f> — » <p + const, which is broken cither 
explicitly or by quantum effects. In the limit of exact 
symmetry, the cj> direction is flat and thus dangerous cor- 
rections to e, -q are controlled by the smallness of the 
symmetry breaking. Moreover, PNGBs like the axion 
are ubiquitous in particle physics: they arise whenever 
an approximate global symmetry is spontaneously bro- 
ken and are plentiful in string compactifications. Axion 
inflation is also phenomenologically desirable since the 
tensor-to-scalar ratio is typically large in such models. 

The first explicit example of axion inflation was the 
natural inflation model [l| in which the shift symmetry 
is broken down to a discrete subgroup <p — > <f> + (2ir)f, 
resulting in a periodic potential 



K P (0)=A 4 [1- cos (<£//)] 



(1) 



with / the axion decay constant. For such potential, 
agreement with observations requires / > M p , which may 
be problematic since it suggests a global symmetry bro- 



ken above the quantum gravity scale, where effective field 
theory is presumably not valid. Moreover, / > M p does 
not seem possible in string theory [§] . More recently, sev- 
eral controlled realizations of axion inflation have been 
studied - including double-axion inflation N-flation 
axion monodromy [f| and axion/4-form mixing Q 

which have / < M p but nevertheless behave effectively 
as large field inflation models ((f) >, M p ). 

In axion inflation models, the inflaton couples to some 
gauge field as ycVF^ 'F^, where F^ u — d^A^ — duA^ and 

p^v = £ m^/3 Fq/3 /2. The scale of this coupling is set by 
the axion decay constant, /; the dimensionless parameter 
a is typically order unity but can be > 1 in multi-field 
or extra-dimensional models Q. It is natural to explore 
the implications of this generic interaction for observ- 
ables. In [7| it was shown that energy dissipation into 
gauge fields can slow the motion of (f>, providing a novel 
new inflationary mechanism that operates at very strong 
coupling. Here, we point out that even in the conven- 
tional slow-roll regime, the coupling (f>FF can have sig- 
nificant impact. The motion of the inflaton amplifies the 
fluctuations of the gauge field, which in turn produce in- 
flaton fluctuations via inverse decay SA + SA — » 5(f). 
When / <, 10~ 2 M p , which is natural for realizations that 
admit an UV completion, we show that the inverse decay 
typically dominates over the usual vacuum fluctuations 
from inflation, and this has dramatic phenomenological 
consequences. Our results are quite general: in the spirit 
of effective field theory, a coupling 4>FF should be in- 
cluded whenever <fi is pseudo-scalar [ill ]. 

Recently, there has been considerable interest in non- 
gaussian effects in the CMB (see the review [l2[ for refer- 
ences). Nongaussianity will be probed to unprecedented 
accuracy with the forthcoming Planck data and may pro- 
vide a valuable tool to discriminate between models. Sev- 
eral constructions are known which can predict an ob- 
servable signature; however, in the minimal cases non- 
gaussianity is small, and obtaining an observable level 
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usually requires either fine-tuning or unconventional field 
theories. Here we point out that the inverse decay con- 
tribution to 66 is highly nongaussian in axion models; 
observational bounds are easily saturated for / not much 
smaller than M p . Thus, the simplest and, perhaps, most 
natural models of inflation can lead to observable non- 
gaussianity. 

II. COSMOLOGICAL PERTURBATIONS 

We consider the theory 

S = -\{d<t>? - V{6) - \F^F, V - ^6F^F^ (2) 

where 6 is the PNGB inflaton, F^ v the field strength of 
the gauge field (for simplicity, a £7(1) gauge field is con- 
sidered; the extension to non-Abelian groups is straight- 
forward), and F^u its dual. The potential V(6) may 
contain a periodic contribution of the form ([T]) due to 
non-perturbative effects and, perhaps, non-periodic con- 
tributions from other effects (such as moduli stabilization 
or wrapped branes). In this Section, we leave V{6) arbi- 
trary, except to suppose that it is sufficiently flat to sup- 
port N e £ 60 e-foldings of inflation. We assume an FRW 
geometry ds 2 = —dt 2 + a(t) 2 dx. 2 — a(r) 2 [— dr 2 + fix 2 ] . 

Working in Coulomb gauge, we decompose A(t, x) into 
circular polarization modes obeying [3] 



dr 2 



k 2 ± 



2fc£ 



A ± (T,k)=0, £ 



2fH 



(3) 



where dot denotes differentiation with respect to t, H = 
a/a, £ = const. We observe that one of the polarizations 
of A experiences a tachyonic instability for kliaH) < 2£. 
The growth of fluctuations is described by |7j 



A+(r,k) = 



1 ( k 



2k \2ZaH 



1/4 



(4) 



in the interval (8£) _1 £ k/(aH) <, 2£ of phase space 
which accounts for most of the power in the produced 
gauge field (we take </> > without loss of generality). 
This interval is nonvanishing only for £ > O(l), which 
we assume in the following. The production is uninter- 
esting at smaller £. 

The unstable growth of A + (r, k) yields an important 
new source of cosmological fluctuations, 66. The pertur- 
bations of the inflaton are described by p], LL3[ 



d 2 . fl V 2 
on — =- 



dt 2 



dt 



66{t,*) = -F^F^ (5) 



where the source term is constructed from (0}. The so- 
lution of (O splits into two parts: the solution of the 
homogeneous equation and the particular solution which 
is due to the source. Schematically, we have 



+ 66 

inv. decay 
homogeneous particular 



(6) 



The quantity of interest is the primordial curvature 
perturbation on uniform density hypersurfaces, ( = 
— ^66. We computed the two-point (C(x)£(y)) an d 
three-point (C (x) C (y) C (z)) correlation functions using 
the formalism of [3, . The two-point function defines 
the power spectrum 



(C(x)C(y)) = / -r 



dk sin [fc|x 
k fc|x — 



y|] 



p c (k) 



We find the result 



P c (k)=V 



1 + 7.5 • iq- 5 v 



(7) 



(8) 



where T 1/2 = 



■, n s is the spectral index, and the 
0.002 Mpc -1 . The two terms in © 



pivot scale is fco 
are the power spectra of the homogeneous and inhomo- 
geneous parts of (O, respectively. There is no "mixed 
term" since the two contributions ([6]) are uncorrelated. 
(The gauge fluctuations that source <5(/>inv. decay, and that 
are amplified according to ([3]) , are not correlated with the 
vacuum inflaton fluctuations.) The power spectrum is 
probed by CMB and Large Scale Structure observations. 
It is found to be nearly scale invariant (n s ~ I, the precise 
value depends on the data set assumed [HI), and have 
amplitude P c {k) = 25 • 10~ 10 [TJ (the so-called COBE 
normalization). When inverse decay fluctuations are sub- 
dominant, we have the standard result V 1 / 2 = 5 • I0~ 5 ; 
however, at large £ the value of V must be modified. 

The three-point correlation function encodes depar- 
tures from gaussianity. Nongaussian effects from inverse 
decays are maximal when all three modes have compara- 
ble wavelength (the equilateral configuration) . The mag- 
nitude of the three-point function is conventionally quan- 
tified using the parameter /jyi, [HI- We find that: 



/•equil 
JNL 



4.4- io 10 :p 



e 6 ^ 
1^ 



(9) 



This result does not include the negligible contribution 
from 66 vac and is accurate as long as |/jvx| !• From 



and © we can see that / 



cquil 
NL 



8400 at large 



cquil 
NL 



< 



eqs. 

f. The current WMAP bounds are -214 < /; 
266 (95% CL), while the Planck satellite, and planned 
missions, will constrain /J^" 1 to O(10) fl6j . 

The results (jl]) and ((9|) only depend on the two dimen- 
sionless combinations £ and V 1 ^ 2 , shown in Figure [TJ The 
solid red curve indicates the parameter values which re- 
produce the COBE normalization of the power spectrum. 
In the region below, and above the dashed black line the 
power spectrum is dominated by 66 v&c and by 66- mv , decay, 
respectively. This second regime is already ruled out by 
the current bound on /^ ul . 

The results ([8} and (|9]) have been obtained by disre- 
garding two backreaction effects of the produced gauge 
quanta. Such quanta are produced at the expense of the 
kinetic energy of 6, so that, if the instability is sufficiently 
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FIG. 1: Values of parameters leading to the observed COBE 
normalization of the power spectrum (red line), and reference 
values for the nongaussianity parameter f^/ = 10, 266, 8000 
along this curve. See the main text for details. 



strong, then it will affect the infiaton dynamics. The re- 
gion of parameter space where this occurs is above the 
black solid line {V 1/2 > 13£ 3/2 e^) shown in Figured] 
We have also disregarded the impact of the energy den- 
sity of the produced quanta on the expansion rate, H. 
This is justified provided e 2 ^/£_ 3 < 2 • 10 4 M 2 /H 2 . This 
constraint is not expressed in terms of £ and V 1 ^ 2 , so 
we have not included it in Figure [1] However, it can be 
studied on a case-by-case basis. 

The gauge quanta also source gravity waves (GW). It 
is customary to normalize the power of GW to that of 
the density perturbations. Proceeding analogously to the 
computation of the density perturbations, we find 



Pgw 



Pc 



.MO 7 



C 



H 

M 1 



2 r 



1 + 4.3 • 1CT 7 



H 2 e 4 < 



mi e 



(10) 



The tensor-to-scalar ratio, r, is an important quantity to 
discriminate between different inflationary models. The 
current observational limit is r <, 0.2 fl4| . and activity is 
underway to probe r > 0.01 [15j . 

III. PREDICTIONS FOR SPECIFIC MODELS 

We now focus our attention on the power-law potential 



V{<t>) = m 4 



(11) 



which subsumes many interesting scenarios. Inflation 
proceeds at large field values <fi > M p and ends when 
(f> ~ M p . For this model, the values of H, <j> and n s 
are uniquely determined by the number of e-foldings of 
observable inflation N e , according to the standard slow 
roll infiaton evolution (e, 77 1) . In the following, we fix 
N e = 60, which is the typical value taken in large field 
models. Once we do so, we are left with the two param- 
eters f/a, and /i. For any given value of f/a, the mass 
scale \i is uniquely determined by fixing the power spec- 
trum ([5]) to the COBE value. We can then plot the other 
observational predictions as a function of f/a only. We 



do so in Figure [2] where we take p — 1, 2 for illustration. 
In both cases, backreaction effects can be neglected. 
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FIG. 2: Observational predictions for the large-field power- 
law inflation model (|ll|l with p = 1,2 and assuming N e = 60. 
The spectral index is n s = 0.975, 0.967 for p — 1,2. At small 
f/a the coupling of 4> to FF is stronger and nongaussianity is 
large. The tensor-to-scalar ratio decreases at strong coupling; 
however, the decrease is important only at values of f/a which 
are ruled out by the current bound on ffi™ . 



Figure [2] shows that large nongaussianity is rather 
generic for large-field axion inflation. The current bound 
is saturated for decay constants / /a < 10~ 2 M p , which is 
natural in a model that admits a UV completion. Cur- 
rent limits on nongaussianity therefore provide an upper 
bound on the strongest couplings of the type <f>FF be- 
tween the infiaton and any gauge field. 

We see also that r decreases at strong coupling. This 
modifies the usual predictions of large field inflation and 
implies, for example, that p — 4 could be made compati- 
ble with observation, at the level of the 2-point function. 

Natural Inflation: The original natural inflation 
model [l[ was based on the potential ((TJ. If we require 
n s £ 0.95, as suggested by recent data [lij], then the 
model requires a large decay constant / > 5M p [I?) • In 
this regime inflation proceeds near the minimum = 
and is indistinguishable from the model (fTTj) with p = 2. 
Large values of / weaken the coupling of <f> to FF, hence 
inverse decay is negligible unless a ^ 200, whereas we 
expect a = 0(1) in the simplest (single-axion) scenario. 
On the other hand, / >, M p may be problematic and it 
seems that a UV completion of axion inflation requires 
/ < M p . We now turn our attention to such scenarios. 

Axion Monodromy: In [5J an explicit, controlled 
realization of axion inflation was obtained from string 
theory. The potential has the form V(4>) — /i 3 + 
A 4 cos(0//) where the linear contribution arises because 
the shift symmetry is broken by wrapping an NS5-brane 
on an appropriate 2-cycle, and the periodic modulation 
is due to nonperturbative effects. The former typically 
dominates [f| |(| so we have the model (ITTj) with p = 1, 
to first approximation. The decay constant is bounded 
[| as 0.06V" 1/2 5 S 1/4 < f/M p < 0.9g s with g s < 1 the 
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string coupling and V 1 the compactification volume 
in string units. From Fig. [5] we see that large nongaus- 
sianity is easily obtained for a — 0(1). Periodic modu- 
lation of V(cj)) can also lead to resonant nongaussianity 
HI for / < 10" 2 Mp and A sufficiently large [al^. 

N-flation: In Q it was noted that the collective mo- 
tion of iV axions fa , each with its own broken shift sym- 
metry, can support inflation when /j < M p , via the as- 
sisted inflation mechanism jl8| . This scenario is quite 
natural in string theory, where generic compactifications 
may contain exponentially large numbers of axions [j-J 0] ■ 
For fa <, fi we can expand the potential near the min- 



r/2. The dynamics of the 



imum to obtain V 
collective field $ = are well-described by a single 

field model of the form (JTTJ with p = 2 0, |||. Sufficient 
inflation requires $ > M p which can be achieved for sub- 
Planckian fa provided y/~N is sufficiently large. 

Typically, the mass basis {fa} is not aligned with the 
interaction basis 0, US] , so that all fa couple to a given 
gauge field as £ int = — Y^i&ifaFF/fi. The coupling of 
the collective field $ to FF is highly mo del- dependent 
but can be parametrized as £j nt = — a c ^F F/f e g, so 
that the result of Fig. [2] apply for the effective coupling. 
A precise calculation depends on the mass rotation and 
the microphysics, but we expect that an observable signal 
is possible for reasonable parameters. 

Double- Axion Inflation: Ref. Q proposed a model 
characterized by two axions 8 and p, with potential 



i=i 



(12) 



which arises from the coupling of the two axions to two 
different gauge groups, j-FjFj, and -^FiFi (up to nu- 
merical coefficients). For fi/gi = /2/S2, one linear com- 
bination of the two axions becomes a flat direction of 
(fT2|) . This relation can be ascribed to a symmetry of the 
theory, and the curvature of the potential along this di- 
rection can be made controllably small if this symmetry 
is only slightly broken. In this case one obtains an effec- 
tive large field inflaton, with a potential of the type ([T}, 
and with an effective axion constant > M p , even if all the 
fi, gi are sub-Planckian. Therefore, this model can lead 



to large production of gauge fields and observable / 



cquil 
NL ■ 



Axion Mixing: Ref. [8j realizes p = 2 via axion/4- 
form mixing. Here / < M p so 3> 1 is possible. 

In summary, we have shown that large nongaussianity 
is possible for many explicit axion inflation models which 
admit a UV completion. Our qualitative results will carry 
over to any inflation model with a PNGB dynamically 
important during inflation, including multi-field models 
such as the roulette [24| or racetrack |25j scenarios. Sim- 
ilar effects may also be possible for higher p-form fields. 
It would be interesting to study the value of a in concrete 
string theory realizations. 
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